Topological entanglement entropy in Gutzwiller projected spin liquids 
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The topological entanglement entropy(TEE) of Gutzwiller projected RVB state is studied with 
Monte Carlo simulation. New tricks are proposed to improve the convergence of TEE, which enable 
us to show that the spin liquid state studied in Ref.[14l] actually does not support Z2 topological 
order, a conclusion that is consistent with the information drawn from the inspection of the topolog- 
ical degeneracy on the same state. We find both a long ranged RVB amplitude and an approximate 
Marshall sign structure are at the origin of the suppression of vison gap in this spin liquid state. On 
the other hand, robust signature of Z2 topological order, i.e., a TEE of In 2, is clearly demonstrated 
for a Gutzwiller projected RVB state on triangular lattice which is evolved from the RVB state 
proposed originally by Anderson[^. We also find that it is the sign, rather than the amplitude of 
the RVB wave function, that dominates the TEE and is responsible for a positive value of TEE, 
which implies that the nonlocal entanglement in the RVB state is mainly encoded in the sign of the 
RVB wave function. Our results indicate that some information that is important for the topologi- 
cal property of a RVB state is missed in the effective field theory description and that a Z2 gauge 
structure in the saddle point action is not enough for the RVB state to exhibit topological order, 
even if its spin correlation is extremely short ranged. 

PACS numbers: 75.10.-b, 73.43.-f, 71.27.+a 



I.INTRODUCTION 

The search for spin liquid phases in quantum antifcr- 
romaget is a hot topic in the field of strongly correlated 
electron systems [ij. A spin liquid state is an exotic state 
of matter in that it can exhibit novel structures and ex- 
citations beyond the Landau- Ginzburg paradigm. The 
topological order and the related fractionalized excitation 
are one of such possibilities. A topological ordered sys- 
tem is unique in that its ground state exhibits topological 
degeneracy on multiply connected manifolds. A topolog- 
ical ordered state is also associated with a nonzero topo- 
logical contribution to the entanglement entropy 0, Si- 

Following the original idea of resonating valence bond 
state, two approaches are widely used for the descrip- 
tion of spin liquid phase. The first one is the effective 
field theory approach based on slave particle representa- 
tion of the spin degree of freedom The second one is 
the variational approach based on Gutzwiller projected 
wave functions j5|]. The Gutzwiller projection procedure, 
which removes the doubly occupied configuration in the 
mean field ground state, amounts to perform gauge av- 
eraging on the saddle point of the effective field theory. 
Thus, if the gauge fluctuation is irrelevant in the long 
wavelength limit, one should expect the Gutzwiller pro- 
jected RVB state to exhibit the same long range physics 
as the effective field theory. In a gapped spin liquid, 
the only nontrivial physics in the long wavelength limit 
is its topological property. However, it is found earlier 
that the topological property predicted by the effective 
field theory does not always agree with that found from 
the corresponding Gutzwiller projected wave function. In 
particular, the Marshall sign rule is found to be closely 
related to the the absence of the topological order in cer- 



tain spin liquid states 

Both the topological degeneracy and the topological 
entanglement entropy(TEE) are used to characterize the 
topological property of a spin liquid state. To detect 
the topological degeneracy, wave functions with differ- 
ent number of trapped visons in the holes of a multiply 
connected manifold are constructed (gI. [lol. Hlj . The over- 
lap of these wave functions, which should decay expo- 
nentially with the linear scale of the system, provides a 
measure of the size of the vison gap(ll|. On the other 
hand, only a single ground state is needed to extract 
the TEEll [H, [iil. The TEE of the Gutzwiller pro- 



jected spin liquid state is first studied in Ref.|l4 Il5 | 
with a Monte Carlo sampling procedure. More specifi- 
cally, they evaluated the TEE of a Gutzwiller projected 
spin liquid state derived from a mean field ansatz with a 
Z2 gauge structure and a very large spin gap, hereafter 
denoted as state RVB-I, and claimed that the obtained 
TEE(0.584± 0.089) is consistent with the expectation for 
a Z2 spin liquid(for which 7 = ln2 w 0.693). However, 
as a result of the poor convergence the true error bar 
is much larger than claimed. A latter calculation with 
higher accuracy by the same group shows that the TEE 
of state RVB-I is far away(0.288 ± 0.107) from the ex- 
pected value of In 2 [HI. 

In this paper, new tricks are proposed to improve the 
convergence of the TEE of the Gutzwiller projected wave 
function in Monte Carlo sampling. With these tricks 
significant reduction in the error bar can be achieved 
and calculation on larger subsystem is possible. We find 
the TEE of the state RVB-I is about 0.2 ± 0.04 for the 
largest subsystem size we have tried, which is signifi- 
cantly smaller than In 2. The TEE is also found to have 
negligible dependence on the size of spin gap. Both of 
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these results indicate that the topological property of 
the state RVB-I is different from the prediction of the 
effective field theory. 

As a crosscheck we have made a test of topological 
degeneracy on the state RVB-I. We find the vison gap of 
this spin liquid state is extremely small, if it does exist at 
all. We find a long ranged RVB amplitude(even though 
the spin correlation is extremely short ranged in the state 
RVB-I) and an approximate Marshall sign structure are 
responsible for the suppression of the vison gap in the 
state RVB-I. We also find that it is the sign of the RVB 
wave function, rather than its amplitude, that dominates 
the TEE and is responsible for a positive value of TEE 
in this state. The reduction of TEE in the state RVB-I 
can thus be understood as a direct consequence of the 
approximate Marshall sign structure in this state. 

With these understandings, we have evaluated the 
TEE of a Gutzwiller projected spin liquid state on trian- 
gular lattice with a short ranged RVB amplitude and 
a frustrated sign. This RVB state can be connected 
smoothly to the short range RVB state first proposed 
by Anderson^, ^IT\ and will be denoted as RVB-II here- 
after. We find the TEE of this state converges steadily to 
the expected value of In 2. Again we find it is the sign of 
the wave function that is responsible for a positive TEE. 
These results indicate that some information that is im- 
portant for the topological property of a RVB state is 
missed in the effective field theory description and that 
the nonlocal entanglement in the RVB state is mainly 
encoded in the sign of the RVB wave function. 

This paper is organized as follows. In the next section, 
we introduce the Gutzwiller projected RVB state and the 
related effective field theory description. We then outline 
the basic steps to evaluate the Renyi entanglement en- 
tropy and the topological entanglement entropy of such 
states. In section III, we introduce new tricks to improve 
the convergence of the TEE. In section IV, we apply these 
tricks to study the topological properties of the spin liq- 
uid state RVB-I. In the same section we also make a test 
of topological degeneracy on the state RVB-I and discuss 
the reason for the suppression of vison gap in this state. 
In section V, we present the result of TEE for RVB-II 
and discuss the relation between the sign structure of 
the RVB wave function and the TEE. In the last section, 
we present some discussions on the relation between the 
effective field theory and the physics of RVB state. 



Here |BCS) is the ground state of the following BCS mean 
field Hamiltonian 

/ Ci+ \ . 

in which ipi — I f I is the Nambu spinor. Uij — 

\ ^ii J 

( Z^*'^ I is a 2 X 2 matrix denoting the hopping 

and pairing order parameters of V'i- On the other hand, 
this mean field Hamiltonian can be viewed as a saddle 
point approximation of a spin liquid state in the effective 
field theory description based on slave Boson representa- 
tion of the spin. The fluctuation of the order parameter 
Vij around their saddle point value takes the form of 
gauge fluctuation in the low energy regime. Depending 
on the structure of C/ij , the gauge fluctuation can be ei- 
ther gapped or gapless in the long wavelength limit. In 
particular, when the gauge symmetry is broken down to 
Z2 in the saddle point, the gauge fluctuation is believed 
to be gapped and the low energy physics of the spin liq- 
uid state is believed to be faithfully represented by the 
saddle point approximation In such a case, the only 
effect of the gauge degree of freedom in the long wave- 
length limit is to induce a Zi topological degeneracy on 
multiply connected manifold. 

From the point of view of the effective field theory, 
the Gutzwiller projection procedure amounts to per- 
form gauge averaging on the saddle point. Thus, when 
the gauge fluctuation is gapped we should expect the 
Gutzwiller projected RVB state to exhibit the same long 
wavelength behavior as the corresponding effective field 
theory. In particular, the topological property predicted 
by both approaches should be the same. 

The topological property of the RVB state can be char- 
acterized by the topological entanglement entropy. To 
define the entanglement entropy on a RVB state, we par- 
tition the system into a subsystem and the corresponding 
environment. The reduced density matrix of the subsys- 
tem is then given by tracing out the degree of freedom in 
the environment 

=Tre I *)(*!, (3) 

and the entanglement entropy of the Von Neumann type 
is defined as 



II. THE GUTZWILLER PROJECTED RVB 
STATE AND ITS TOPOLOGICAL 
ENTANGLEMENT ENTROPY 

The Gutzwiller projected BCS mean field ground state 
is widely used to describe a spin liquid state. It has the 
following general form 

14-) =Pg|BCS). (1) 



S = -Trps Inps. (4) 

The entanglement entropy in general satisfies the area 
law S = aL — 7, in which L is the length of the boundary 
separating the subsystem from the environment. 7 is 
the universal topological contribution and is nonzero in 
a topological ordered state. 

It is in general quite hard to evaluate the Von Neu- 
mann entropy from the definition. However, it is proved 
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that the TEE can also be extracted from the more gen- 
eral Renyi entropy [l8|. The Renyi entropy of order n is 
defined as 



Sn 



1 



1 



InTrp^. 



(5) 



The Von Neumann entropy is related to the Renyi en- 
tropy hy S = lim„_).i Sn- It is proved that the topolog- 
ical contribution to the Renyi entanglement entropy is 
independent of the value of n[l8l|. 

In the following, we will use S2 to extract the TEE. 
To evaluate S2, we use the replica method and introduce 
an identical copy of the system 14 , [l^ . The replica is 
partitioned in exact the same manner as the system into 
a subsystem and an environment. Denoting the state 
vector of the system and its replica as a whole as |^'(g)\I'), 
it can shown that 



Trpl = (* ® ^'|SWAP|* ® *) 



(6) 



Here SWAP denotes the operation of exchanging all the 
degree of freedoms within the subsystem between the sys- 
tem and its replica[l3. To evaluate Trp^, we introduce an 
orthonormal basis {si\ in both the system and its replica 
and denote the wave function in this basis as ip{si). Then 
we have 



Trp^ 



2 i^{si)ip{s-, 



(7) 



here Si and Si' are the spin configurations of the system 
and its replica before the swap operation, s; and Si' are 
the swapped spin configurations in the system and its 
replica. In Ref.[l4|, a sign trick is proposed to evaluate 
the above average. For this purpose we rewrite the above 
average as follows 



in which 

(SWAPamp 



(SWAP) = (SWAPanip) X (SWAP,ig„), 



(8) 



(SWAPs 



(9) 



Here p{s„s,') = |V'(Si)P|V'(si')P> p(si,Si' 
|V'(si)V'(si')V'(sj)V'(si')l- (t>(si,Si') is the phase of 



V>(si)Vi(s^/) 



For both the spin liquid state RVB-I and 



'/'(si)V'(Si') ' 

RVB-II studied in this work, it can be shown that the 
wave function '4'{si) is real up to a global phase as a re- 
sult of the time reversal symmetry. Thus e**^*^**'**'-* — ±1. 

In general, we expect both (SWAPsign) and 
(SWAPamp) to decrease exponentially with the sub- 
system size, as is required by the area law for the 
entanglement entropy. Following Ref . (l6| . we define 



^sign 



and 



S'™'' as the contribution of the phase and 



the amplitude of the wave function to the entanglement 
entropy. They are given by 



oamp 

•^2 



ln(SWAP,ig„) 
= - ln(SWAPa,„p) 



(10) 




FIG. 1: An illustration of the partition of the degree of free- 
doms on the square lattice used to extract the TEE from 
Eq.(ll). 

In principle, the TEE can be extracted from the defin- 
ing equation 5*2 = oiL — 7 by an appropriate extrapola- 
tion procedure. However, the result thus obtained suffers 
from the ambiguity in the definition of L on a lattice. To 
remove such ambiguity, a bulk and boundary cancelation 
procedure is widely used. As in Rcf.jlJ]) we divide the 
subsystem into three parts. A, B and C. An illustration 
of the partition on the square lattice is shown in Fig.l. 
Then the topological contribution to the entanglement 
entropy can be extracted by subtracting out both the 
the bulk and boundary contributions from the following 
combination 



QlJ qAU qAU qUIJ I c 
"-"2 ^ "^2 ^ '-'2 ~ "-"2 "T '^2 



ABC 



(11) 



III. THE TRICK OF RE- WEIGHTING AND THE 
TRICK OF RATIO ESTIMATOR 



The evaluation of TEE involves the computation of 
= the sign contribution 82^''^ and amplitude contribution 
82"^^ for each of the seven subsystems in Eq.(ll). For 



a general RVB state, both iSj'^" and S'™'' follows the 
area law. However, the logarithmic correction to the area 
law in a gapless spin liquid state is found to be mainly 
contributed by S'2'*^"[I^. 

As we will see below, for the state RVB-I, 6*2 is mainly 
contributed by S'™^ as a result of an approximate Mar- 
shall sign rule. For large subsystem size, the conver- 
gence of 5*2™^ becomes very slow. This is expected since 
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the spin configurations before and after the swap opera- 
tion become increasingly remote from each other in the 
Hilbert space with the increase of the subsystem size. As 
a result, the fluctuation in the quantity to be averaged, 
R = ^(s'IvjI^'o ' g'^O'^s rapidly with the subsystem size. 
However, the averaged value of such a ratio should de- 
crease exponentially with the length of the boundary, as 
is required by the area law. A Monte Carlo sampling of 
such a wildly fluctuating quantity is thus very hard and 
this explains why Ref.[14] fails to reach the true conver- 
gence in the value of TEE. 

In the following, we propose two tricks to improve the 
convergence of the 82"^^- The first trick is a simple re- 
weighting procedure. Rather than sampling the weight 
W{si,Si') — p{si,Sii), we sample the combined weight 
W'{si,Si') — p{si,Si') + p{si,Si'). It is shown in Ref.[20[ 
that such a combined sampling can reduce the statistical 
fluctuation. With such a combined weight, the quantity 
to be averaged becomes 



(SWAPamp) 



j:s^„ w'{s.,s,)j^ ^ ^^^^ 



Now the quantity to be averaged becomes 
whose fluctuation is the square root of R. 

For even larger subsystem size, the above re- weighting 
trick will also fail as the fluctuation in R grows expo- 
nentially with the subsystem size. A trick to cope with 
such an exponentially growing fluctuation is the ratio 
estimator [19|, [21 1. In such a trick, we express the quan- 
tity to be evaluated as the product of a series of ratios 
so that the evaluation of each ratio only suffers from a 
much smaller fluctuation. There are different ways to 
implement this trick. One convenient way is to rewrite 
(SWAP amp) in the following form 



(SWAP 



A 



amp/ 



n 



y 



(13) 

in which A = \ip{si)ip{si')tp{si)il:{si>)\ , B = 
|^(si)PlV'(si')P- ^ [0; 1] is a series of powers satis- 
fying ri < Ti^i, ri = and r„i^i = 1. If we define 
Wi{si,Si') — yl^^'"'+i_B'''+i , then it can be shown that 



IV. THE TOPOLOGICAL PROPERTIES OF THE 
STATE RVB-I 

In Ref.^, the TEE of the spin liquid state with the 
following Zi mean field ansatz is studied. 



= ry(ri-f 



(15) 



in which r^'^''^ are the three Pauli matrices, 77 and A 
are two real parameters. A is determined self consis- 
tently from the mean field condition of {ciiCi-\) = 0. The 
Gutzwiller projected RVB state derived from this ansatz 
describes a symmetric spin liquid. The RVB amplitude 
is given by 



a(Ri — Rj 



ik-(Ri-Rj 



k 



(16) 



in which = — 2(cos(ka;) + cos(k3;)), Ak = A -I- 
47ycos(ka;) cos(ky) -I- 4i77sin(ka;) sin(ky). At the effective 
field theory level, this spin liquid phase is predicted to 
possess Z2 topological order since the gauge symmetry 
is broken down to Z2 at the saddle point and that spin 
excitation spectrum is fully gapped. We thus expect to 
see a TEE of In 2 in the state RVB-I. 




FIG. 2: The Renyi entanglement entropy from the amplitude 
and the sign contribution in the state RVB-I. The subsystem 
is a 1/ X L square region and the total system size is 16 x 16. 



When Ti+i — Ti are chosen to be sufficiently small, each 
term in the above product can be evaluated easily. We 
note a similar trick is recently proposed to compute the 
TEE of quantum dimer model at the R-K point by the 
present authors 21 1. 



In our calculation of the entanglement entropy, we find 
the sign contribution 53'^" for this state is much smaller 
than the amplitude contribution S*™^, as is illustrated 
in Fig. 2. The origin for such an unusual behavior is an 
approximate Marshall sign structure in this RVB state. 
Although the ansatz Eq.(15) is frustrated and the Mar- 
shall sign rule is in a strict sense violated, we find the 
extent to which the Marshall sign rule is violated is al- 
ways limited. To illustrate this point, we plot in Fig. 3 the 
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average of the Marshall sign in the state RVB-I. The av- 



eraged Marshall sign is defined as ((— 1) '^)[22|, in which 



is the number of down spins in sublattice A. 




FIG. 3: The averaged Marshall sign in the spin liquid state 
RVB-I as a function of the gap parameter rj. 

As a result of such a sign rule, 52*^" converges rapidly. 
However, the evaluation of is very challenging. 

In the evaluation of S'™^, we have used both the re- 
weighting and the ratio estimator trick introduced in the 
last section. In Fig. 4a, we show the result of TEE for 
various subsystem sizes. The TEE is seen to be much 
smaller than the value of In 2 predicted by the Z2 effec- 
tive field theory. To test if the deviation from the effective 
field theory prediction is caused by the finite size effect in 
spin correlation, we have studied the dependence of TEE 
on the size of the spin gap, the result of which is shown 
in Fig. 4b. The TEE is found to be insensitive to the size 
of spin gap. Thus it is unlikely that the deviation from 
effective field theory prediction is caused by a finite size 
effect in the spin correlation. 

As a further check of the above result, we perform a 
test of topological degeneracy on the state RVB-I. A spin 
liquid state with Z2 topological order should exhibit four- 
fold topological degeneracy on a torus. The four degener- 
ate states differ from each other in the number of trapped 
visons in the two holes of the torus. To detect such de- 
generacy, we evaluate the overlap of the wave functions 
with different number of inserted visons in the holes of a 
finite torus. In a topological ordered state, such overlaps 
should vanish in the thermodynamic limit. The decay 
rate of the overlap provides us an estimate of the size of 
the vison gap[llj|. A nonzero vison gap is the defining 
property of a spin liquid state with Z2 topological order. 
It should be noted that the vison gap has nothing to do 
with the spin gap. In fact, a spin liquid state can have a 
zero vison gap even when its spin correlation is extremely 
short ranged. A well known example of this kind is the 
short range RVB state defined on the square lattice, for 
which the vison gap vanishes as a result of the bipartite 





FIG. 4: (a)The TEE of the state RVB-I as a function of the 
subsystem size. Here La is the linear scale of region A in 
Fig.l. The total system size is fixed at 16 x 16. (b)The TEE 
as a function of the gap parameter rj. The total system size 
is fixed at 12 x 12 and the subsystem size is fixed at La ~ 3. 



nature of the lattice[23|. 

For the Gutzwiller projected RVB state, inserting a 
Z2 vison in a hole of torus is equivalent to changing the 
boundary condition of the RVB amplitude around that 
hole from periodic to anti-periodic or vice versa. For 
convenience's sake, we choose to evaluate the overlap be- 
tween the RVB state with periodic-antiperiodic boundary 
condition and that with antiperiodic-periodic boundary 
condition. In the calculation we have set 77 = 1.5 as in 



Ref. IJ], which corresponds to a spin correlation length 
as short as 1.3 lattice constant. The result of the overlap 
is shown in Fig. 5. It fluctuates with the system size and 
decay only very slowly if compared with the decay rate 
of the spin correlation function. This would imply an ex- 
tremely small vison gap, if it does exist at all. We thus 
find the result of the topological degeneracy is consistent 
with result of TEE and both indicate that the topological 
order in the state RVB-I is very weak, if it exist at all. 
The result in Fig. 5 also indicates that to get a converged 
result for TEE, one should use subsystem with size much 
larger than 20 x 20, which is of course unrealistic. 

Thus for the spin liquid state RVB-I, the effective field 
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FIG. 5: The overlap of the wave functions with different num- 
ber of inserted visons in the holes of a torus, (a) in linear scale 
(b) in logarithmic scale. The decay rate of the overlap is pro- 
portional to the inverse of the vison gap [Til]. 



theory and the Gutzwiller projected wave function pre- 
dict totally different topological properties. To under- 
stand the origin of such a discrepancy, we note that it 
is previously found that the Marshall sign rule in the 
RVB wave function can result in the suppression of the 
vison gap and can lead to the absence of topological 
degeneracy [^-Q. For similar reasons, the approximate 
Marshall sign structure in the state RVB-I will definitely 
enhance the overlap integral and thus reduce the vison 
gap. 

There is another reason for the suppression of vison gap 
in the state RVB-I. As is well known, the Z2 topological 
classification of a RVB state is only meaningful when the 
RVB amplitude is short ranged. However, the RVB am- 
plitude for the state RVB-I is long ranged, even though 
its spin correlation is extremely short ranged. This is 



because the integrand 



in Eq.(16) is singu- 



lar when |Ak| — and < 0. From the self-consistent 
equation (c^Ci-t-) = 0, we have 



It can be easily shown that to satisfy this equation Ak 
must cross zero between the F — (0, 0) and the M = (0, tt) 
point in the Brillouin zone. Thus the gap node must 
occur inside the Fermi surface and RVB amplitude of 
the state RVB-I can not be short ranged, although the 
spin excitation spectrum is fully gapped. Similar situa- 
tion also occurs in a p -|- ip'-wave BCS superconductor, 
in which the gap node at the F point forces the wave 
function of the Copper pair to have a long range tail, al- 
though the Fermi surface is fully gapped by the complex 
pairing |24j. 



V. THE TOPOLOGICAL PROPERTIES OF THE 
STATE RVB-II 

From the results of the last section, we see the 
Gutzwiller projected RVB state may have different topo- 
logical property from that predicted by the corresponding 
effective field theory, even if the gauge symmetry is bro- 
ken to Z2 and the spin excitation is fully gapped. More 
specifically, we find both the Marshall sign rule and a 
long ranged RVB amplitude can act to suppress the vi- 
son gap. Following this line of reasoning, we expect that 
the Z2 topological order predicted by the effective field 
theory should be robust in a Gutzwiller projected RVB 
state when it has a short ranged RVB amplitude and a 
frustrated sign. 

To illustrate this point, we consider a spin liquid state 
defined on the triangular lattice. It has the following 
mean field ansatz 



0. 



(17) 



(18) 



Here ai.2,3 are the three elementary translation vectors 
on a triangular lattice, ii is the coordinate of a lattice 
site in the oi direction. The sign of the pairing terms are 
such that around each elementary rhombic plaquette of 
the triangular lattice the product of the pairing term is 
negative. This ansatz is illustrated in Fig. 6. 

This state is first studied in Ref.|5|]. When ^ — > — cxj, 
the Gutzwiller projected RVB state derived from the 
above ansatz evolves into the short ran ge R VB state pro- 
posed by Anderson on triangular lattice 17| . On the other 
hand, when /i 0, it becomes a rather good variational 
ground state of the Heisenberg model on triangular lat- 
tice and possesses Dirac type spinon excitation at low 
energy. The RVB state for general value of /i can be 
viewed as an interpolation between these two limits. 

For nonzero value of /i, the gauge symmetry of the 
ansatz is broken to Z2 and the spin excitation is fully 
gapped. At the same time, the RVB amplitude is short 
ranged and the sign of the wave function is frustrated. So 



FIG. 6: An illustration of the triangular lattice and the mean 
field ansatz Eq.(18). The pairing term is positive on the red 
bonds and negative on the black bonds. 

we have every reason to expect a robust Z2 topological 
order in this state. 




L 



FIG. 7: The contribution from the sign and amplitude of the 
wave function to the Renyi entanglement entropy in the state 
RVB-II. The subsystem is now a L x L rhombic region on the 
triangular lattice and the total system size is fixed at 20 x 20. 

In the calculation we set /i — —30 to produce a rather 
short ranged RVB amplitude. Unlike the state RVB-I, 
the state RVB-II has a fully frustrated sign. We find 
the contribution from the sign now makes up a much 
larger portion in the total entanglement entropy. This is 
illustrated in Fig. 7. The result for TEE is shown in Fig. 8. 
The TEE is seen to converge steadily to the expected 
value of In 2 for a spin liquid with Z2 topological order. 

In a previous studv|16|. it is found that in a critical 
RVB state the logarithmic correction to the area law is 
mainly contributed by the sign of the wave function. To 
see how the sign of the wave function would contribute to 
TEE, we have plotted in Fig. 9 the contribution to TEE 



FIG. 8: The TEE of the state RVB-II as a function of the 
subsystem size. Here La is the linear scale of a rhombic region 
on triangular lattice. The total system size is fixed at 20 x 20. 

from the sign and the amplitude of the wave function 
separately. A key observation is that for both the state 
RVB-I and RVB-II, the contribution from the sign and 
that from the amplitude always have opposite signs. It is 
always the sign, rather than the amplitude of wave func- 
tion that generate a positive contribute to TEE. It is also 
interesting to note that although 6*2™^ > 6*2'^" in both 
RVB-I and RVB-II, it is the contribution from the sign 
that dominate the TEE. This indicate that the topolog- 
ical property of the RVB state is mainly determined by 
its sign structure. 

VI. DISCUSSIONS 

In this paper, we have evaluated the TEE for two 
Gutzwiller projected Fermionic RVB states. Both states 
are derived from mean field ansatz with a Z2 gauge struc- 
ture and a full gap to spin excitations. Following effec- 
tive field theory arguments, both states should exhibit Z2 
topological order. However, from the result of TEE we 
find the two states have different topological properties. 
With new tricks to improve the convergence of TEE, we 
are able to show that the spin liquid state RVB-I studied 
in Ref. does not support Z2 topological order, while 
the spin liquid state RVB-II proposed in Ref.Q exhibits 
robust Z2 topological order. 

These results indicate that something beyond the effec- 
tive field theory is important to the topological property 
of a RVB state. We find both a long ranged RVB ampli- 
tude and an unfrustrated sign in the RVB wave function 
can act to suppress the vison gap and thus the topologi- 
cal order in a RVB state. It is important to note that for 
the Fermionic RVB state, the RVB amplitude need not 
be short ranged even if the system has a large spin gap 
and that its spin correlation is extremely short ranged. 
The gap node deep inside the Fermi surface can have 
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FIG. 9: The contribution of tiie sign and the amphtude of the 
wave function to the TEE in the state RVB-I(a) and RVB- 
11(b). The total system size is fixed at 16 x 16 and in (a) and 
is fixed at 20 x 20 in (b). 



structure and a fully gapped spin excitation spectrum 
is enough to guarantee the stability of the saddle point 
against gauge fluctuation up to Gaussian level. Since 
Gutzwiller projection just amounts to gauge averaging 
on the saddle point, we expect the Gutzwiller projected 
RVB state to have the same topological property as the 
saddle point action of effective field theory, if the fluctua- 
tion effect beyond the Gaussian level can be neglected in 
the long wavelength limit. Our results show that this is 
not always true. Thus, certain singular fluctuation mode 
beyond the Gaussian level is important for the topolog- 
ical property of RVB state. It is an interesting prob- 
lem to study how such singular gauge mode are related 
to the emergence of the Marshall sign rule structure of 
the Gutzwiller projected wave function and to the long 
ranged nature of the RVB amplitude. 

In conclusion, from our results it is clear that the effec- 
tive field theory at the Gaussian level misses two pieces 
of information that is important for the topological prop- 
erty of a RVB state. The first one is the sign structure of 
the Gutzwiller projected wave function. The second one 
is the range of the RVB amplitude of the RVB state. The 
effective field theory can predict the topological property 
of a RVB state reliably only when both of these informa- 
tion are accounted for. 

This work is supported by NSFC Grant No. 10774187, 
No. 11034012 and National Basic Research Program of 
China No. 2010CB923004. 



profound effect on the long wavelength behavior of the 
system. This is markedly different from a Bosonic RVB 
state, for which a full spin gap always corresponds to a 
short ranged RVB amplitude. 

We have also tested topological degeneracy on the state 
RVB-I. We find the result implies a very small or even 
zero vison gap for the state RVB-I, which is consistent 
with the conclusion we get from the result of TEE. To our 
knowledge, this is for the first time that the topological 
property of a RVB state is studied from both the TEE 
and the topological degeneracy perspective with consis- 
tent conclusions reached. We also find that it is the sign, 
rather than the amplitude of the wave function that dom- 
inates the TEE and is responsible for a positive value of 
TEE in the Gutzwiller projected RVB state. Since a 
positive TEE is the signature of nonlocal entanglement 
in the system, we find the nonlocal entanglement in a 
RVB state is mainly encoded in its sign structure. This 
is consistent with previous studies in which it is found 
that a Marshall sign rule is at the origin for the absence 
of topological order in certain spin liquids d, 0, l^- 

From the effective field theory point of view, a Z2 gauge 
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